Probing dark energy using baryonic oscillations in the galaxy 
power spectrum as a cosmological ruler 
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ABSTRACT 

We show that the baryonic oscillations expected in the galaxy power spectrum may be used as 
a "standard cosmological ruler" to facilitate accurate measurement of the cosmological equation 
of state. Our approach involves a straight-forward measurement of the oscillation "wavelength" 
in Fourier space, which is fixed by fundamental linear physics in the early Universe and hence is 
highly model- independent. We quantify the ability of future large-scale galaxy redshift surveys 
with mean redshifts z ~ 1 and z « 3 to delineate the baryonic peaks in the power spectrum, 
and derive corresponding constraints on the parameter w describing the equation of state of the 
dark energy. For example, a survey of three times the Sloan volume at z « 1 can produce a 
measurement with accuracy Aw w 0.1. We suggest that this method of measuring the dark 
energy powerfully complements other probes such as Type la supernovae, and suffers from a 
different (and arguably less serious) set of systematic uncertainties. 

Subject headings: cosmological parameters — large-scale structure of universe — surveys 



1. Introduction 

Measurement of anisotropies in the Cosmic Mi- 
crowave Background (CMB) radiation have shown 
the Universe to be (very close to) spatially flat. 
However, matter only makes up a third of the crit- 
ical density - the dominant contribution to the 
energy density of the Universe appears to exist in 
an unclustered form called "dark energy" . Fur- 
thermore, observations of distant supernovae have 
indicated that the expansion of the Universe has 
entered a phase of acceleration. This is explica- 
ble if dark energy has recently come to dominate 
the dynamics of the Universe and exerts a large 
negative pressure. This situation occurs naturally 
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if Einstein's equations contain a cosmological con- 
stant, but the observed magnitude of the vacuum 
energy is wildly inconsistent with the current pre- 
dictions of quantum physics. This has motivated 
the consideration of a more general dark energy 
equation of state, P — w p (Turner & White 1997). 
An accelerating universe is produced if w < — 1/3 
and vacuum energy is described by w = — 1. If the 
dark energy is the product of new physics, such as 
quintessence (e.g. Ratra & Peebles 1988), then 
the resulting equation of state varies with redshift 
and must be treated as a more general w(z) (e.g. 
Linder & Huterer 2003). 

Securing more accurate measurements of the 
dark energy component is of prime cosmological 
importance. A number of methods have been in- 
vestigated, most notably the use of Type la su- 
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pernovae as "standard candles" (e.g. Weller & 
Albrecht 2002). Other possible probes include 
counts of galaxies (Newman & Davis 2000) and of 
clusters (Haiman et al. 2001), weak gravitational 
lensing (Cooray & Huterer 1999), the Alcock- 
Paczynski test applied to small-scale galaxy cor- 
relations (Ballingcr et al. 1996) and use of the 
CMB (e.g. Douspis et al. 2003). In this paper 
we examine a probe of dark energy that has not 
received much attention in the literature (e.g. La- 
hav 2002), but which we believe has many advan- 
tages: baryonic oscillations in the galaxy power 
spectrum. 

The CMB power spectrum contains acoustic 
peaks. Coupling between baryons and photons at 
recombination imprints these "wiggles" into the 
matter power spectrum on a scale corresponding 
to the sound horizon in the early Universe (Peebles 
& Yu 1970, Eisenstein & Hu 1998). The positions 
of the peaks and troughs in Fourier space are cal- 
culable from straight-forward linear physics and 
act like a "standard cosmological ruler" (Eisen- 
stein, Hu & Tegmark 1998). Therefore a power 
spectrum analysis of a galaxy redshift survey 
containing acoustic oscillations can be used to 
measure the cosmological parameters (Eisenstein 
2002): the conversion of the redshift data into 
real space requires values of the parameters to be 
assumed; an incorrect choice leads to a distor- 
tion of the power spectrum and the appearance of 
the acoustic peaks in the wrong places. Standard 
ruler techniques for deducing the cosmic parame- 
ters have been proposed many times before. In the 
Alcock-Paczynski test (Alcock & Paczynski 1979) 
a spherical system is distorted into an ellipsoid in 
the wrong cosmological model. Other examples 
include the detection of a characteristic cluster- 
ing scale within the 2dF quasar survey (Roukema, 
Mamon & Bajtlik 2002) and the use of the angular 
power spectrum of dark matter haloes (Cooray et 
al. 2002). Our study is focussed on the ability of 
the standard ruler provided by the baryonic peaks 
to measure the dark energy. 

Mapping the acoustic peaks in the galaxy power 
spectrum at high precision, matching those al- 
ready measured in the CMB power spectrum, 
would provide a spectacular confirmation of the 
standard cosmological model in which mass over- 
densities grow from the seeds of CMB fluctua- 
tions. The sharp features of the baryonic oscil- 



lations represent a powerful and precise observa- 
tional test of the current cosmological paradigm 
("A-CDM"). Accurate delineation of the baryonic 
oscillations lies beyond the current state-of-the-art 
galaxy redshift surveys, the 2dF Galaxy Redshift 
Survey (Colless et al. 2001) and Sloan Digital Sky 
Survey (York et al. 2000). Moreover, it is advan- 
tageous to place the survey volume at higher red- 
shift (Eisenstein 2002) , where the linear regime of 
galaxy clustering extends to smaller physical scales 
(unveiling acoustic oscillations to higher spatial 
frequencies) . 

In this initial study we present a simple "proof- 
of-concept" that acoustic oscillations measured at 
high redshift may be used to accurately deter- 
mine the dark energy parameter w, restricting our 
attention to cosmological models with constant 
(redshift-independent) w. We perform simulations 
to quantify the scale of redshift survey required to 
recover the wiggles and determine the resulting 
constraints on w. We suggest that measurement 
of the dark energy should be a prime motivation 
for the next-generation galaxy redshift survey. 

2. Measuring baryonic oscillations 

In this section we investigate the scale of red- 
shift survey required to measure the acoustic oscil- 
lations in the galaxy power spectrum with a spec- 
ified accuracy. 

2.1. Assumptions and approximations 

We modelled the baryonic oscillations using the 
transfer function fitting formulae of Eisenstein & 
Hu (1998). The matter power spectrum P(k) 
is deduced from the transfer function T(k) us- 
ing P(k) = Ak n T 2 (k); we took the primordial 
spectral slope to be n — 1, as (approximately) 
suggested by inflationary models, and fixed the 
normalization A such that erg — 1. The transfer 
function model is specified by assigning values to 
the matter density fi m , the baryon density Clb and 
Hubble's constant h = H /(W0 km s" 1 Mpc" 1 ). 
We assumed that it does not depend on the dark 
energy equation of state: the energy density of 
dark energy, p Xl relative to that of matter, p m , 
scales with redshift as p x I Pm oc (1 ~t~ z^ w and is 
a negligible contributor to physics before recom- 
bination (as w < — 1/3). Of course more com- 
plex dark energy models with varying w may lead 



2 



to significant effects at the epoch of recombina- 
tion; we do not consider these here, but point out 
that the CMB is well described by simple matter- 
dominated models (Spergel et al. 2003). In this 
section we assume Q m — 0.3, Slb/Cl m — 0.15, 
h = 0.7 and w = — 1 (i.e. a cosmological con- 
stant). Our assumed value for fit, is consistent 
with the latest CMB results (Spergel et al. 2003), 
constraints from Big Bang nucleosynthesis theory 
(e.g. O'Meara et al. 2001), and the shape of 
the galaxy power spectrum measured by the 2dF 
Galaxy Redshift Survey (Percival et al. 2001). We 
assume throughout this study that the Universe is 
flat, fi tot = 1. 

The baryonic oscillations can be conveniently 
displayed by dividing the model P(k) by a smooth 
reference spectrum containing no wiggles (lower 
panel of Figure 1; the reference spectrum was ob- 
tained from the fitting formulae of Eisenstein & Hu 
(1998) together with the transfer function). The 
result is well-approximated by a slowly-decaying 
sinusoidal function whose peaks are harmonics of 
the sound horizon scale: approximate formulae for 
their positions in Fourier space are given by Eisen- 
stein & Hu (1998). The amplitude of the wiggles 
increases with the baryon density fif,. The key 
quantity we wish to observe is the "wavelength" 
of the acoustic oscillations in fc-space, which we 
denote as kA, and which is related to the sound 
horizon at recombination by kA = 2ir/s. 

If dark energy is neglected at high redshift, the 
co-moving sound horizon size at last scattering, s, 
is given by (Hu & Sugiyama 1995, Cornish 2001): 

where a r , a cq are the values of the scale fac- 
tor a = 1/(1 + z) at recombination and matter- 
radiation equality respectively. c s is the speed of 
sound and is approximately over the interval 
of integration. The value of s is of order 100 ft -1 
Mpc. 

Thus the theoretical value of Ua is set by fun- 
damental CMB physics and depends strongly on 
fi m , weakly on fij, and negligibly on dark energy. 
Ha is our "standard measuring rod" . In a redshift 
survey at intermediate redshift z, the apparent size 
of the oscillation wavescale Ua will depend on the 
cosmological geometry, now including the effects 
of dark energy. As we will see in Section 3, the 
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Fig. 1. — The model power spectrum of Eisen- 
stein & Hu (1998) for cosmological parameters 
fi m = 0.3, f2b/f2 m = 0.15, ft = 0.7. For this choice 
of parameters, the sound horizon s m 105 ft -1 Mpc 
and the wiggle wavescale Iza ~ 0.0601ft- Mpc -1 . 
In the upper panel we divide the model P(k) by 
the corresponding zero-baryon model, fi m = 0.3, 
fib = 0, ft = 0.7. Replacing cold dark matter 
with baryons produces acoustic oscillations and an 
overall suppression of power. In the lower panel we 
divide the model P(k) by a smooth fit to the over- 
all shape of the spectrum. The arrows indicate 
the approximate position of the linear/non- linear 
transition at different redshifts, estimated in the 
following way. At z = 0, we conservatively de- 
fined the linear regime by k < k D \ = 0.1 ft Mpc -1 . 
From the model P(k) we computed the variance 
of mass fluctuations <J 2 (R) inside a sphere of ra- 
dius R, where R = ir/2k n \ (i.e. half a fluctuation 
wavelength, or whole wavecrest, coincides with the 
diameter 2R). For k n \ = 0.1 h Mpc -1 we found 
<7 2 (7r/2fc n i) = 0.35, and then applied this criterion 
to fix the linear/non- linear transition at other red- 
shifts. At higher z, the amplitude of P(k) is re- 
duced by the growth factor, P(k) — > P(k) Di(z) 2 . 
At z = 1, for example, a 2 (n/2k) = 0.35 for k = 
0.19 ft Mpc -1 = k nl (z = 1). We fixed the overall 
amplitude of P(k) such that cr 2 (8ft -1 Mpc) = 1. 
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effects of assuming an incorrect world model (e.g. 
incorrect w) would include a distortion of the mea- 
sured value of kA ■ (There will be different effects 
parallel and perpendicular to the line-of-sight, but 
these are averaged out in our approach) . To zeroth 
order, the precision with which we can empirically 
measure Ua tells us how accurately we can mea- 
sure the geometrical distance to the redshift z and 
hence how accurately we can measure the equation 
of state. 

We have assumed that the power spectra of 
galaxies and of matter are related by a linear bias 
factor b (where P ga i = b 2 P ma t)- Provided our sim- 
ulated surveys are dominated by sample variance 
(see Section 2.2), the value of b is unimportant: 
the fractional error bars opjP (i.e. the appear- 
ance of Figures 1 and 2) do not change if P{k) is 
scaled by a constant factor (see equation 2). How- 
ever, increasing the value of b (hence increasing 
Pgai) reduces the fractional error due to shot noise, 
and hence reduces the number of objects required 
to render shot noise negligible, by a factor b 2 . The 
assumption of linear bias is simplistic and we plan 
to include more complicated (redshift-dependent) 
biassing schemes in a future study. 

We scaled the power spectrum to redshift z 
using the appropriate growth factor, P(k, z) = 
P(k, 0) Di(z) 2 . Throughout this paper we use the 
approximation of Carroll, Press & Turner (1992) 
for Di(z). This formula is only valid for w = — 1, 
but is a satisfactory estimate for models with 
w ^ — 1, given that we only treat small depar- 
tures from w = — 1 and that we divide out the 
overall shape of P(k) in the analysis procedure. 

These model power spectra are only valid in the 
linear regime of structure formation; on smaller 
scales the non-linear growth of structure washes 
out the baryonic oscillations. At redshift zero, the 
linear regime only extends to k < 0.1 h Mpc -1 
(see Meiksin, White & Peacock 1999, Percival et 
al. 2001) and perhaps only the first major peak 
shown in Figure 1 is preserved (i.e. the peak at 
k rts 0.075 h Mpc -1 ; there is a small peak predicted 
at k w 0.025 h Mpc -1 which is undetectable due 
to cosmic variance). At higher redshifts [z > 1) 
the linear regime extends to much smaller scales, 
unveiling more acoustic peaks. Moreover, a high- 
redshift study is necessary to probe dark energy 
effectively: the conversion of redshifts to spatial 
scales only depends on the Hubble constant H 



at low z and cannot distinguish between different 
values of w. 

A galaxy power spectrum deduced from a red- 
shift survey is subject to redshift-space distortions 
owing to the peculiar velocities galaxies possess on 
top of the bulk Hubble flow. There are two effects 
(Kaiser 1987): the coherent infall of galaxies into 
concentrations of mass (which boosts power on all 
scales by a constant factor) and the incoherent ve- 
locities of galaxies in the central regions of clusters 
(which damps power on small scales). We make no 
attempt to simulate these effects. Rcdshift-space 
distortions will only produce smooth changes in 
the shape of the angle-averaged P(k), not sharp 
features such as acoustic peaks, and we note again 
that the fractional error bars crp / P do not depend 
on the absolute value of P(k) (equation 2). In 
fact, we explicitly divide out the shape using the 
smooth reference spectrum. Moreover, incoherent 
velocities do not have an important effect on the 
large scales probed by acoustic oscillations. These 
velocities can be modelled by a ID pairwise disper- 
sion <7p, and become important on scales k where 
k <T p > H n . At redshift zero, <j v w 300 km s -1 and 
power is hence damped on scales fc > 0.3 /i Mpc -1 . 
Furthermore, simulations indicate that the value 
of <jp drops by a factor ~ 2 between z — and 
z = 2 (Magira, Jing & Suto 2000). 

Some models of galaxy clustering (Peacock & 
Smith 2000, Seljak 2000) predict that the baryonic 
oscillations in the power spectrum are diluted by 
a wiggle-free halo contribution (which amounts to 
a non-local recipe for bias, in which the probabil- 
ity of finding a galaxy is not a simple function of 
the local mass density). In our initial treatment 
we assume a standard linear bias factor. However, 
this should not modify our results to first order: at 
z = 0, the halo contribution to the matter power 
spectrum becomes equal to the linear contribu- 
tion at k w 0.4 h Mpc -1 (Peacock & Smith 2000, 
Figure 2), whereas most of our high-redshift con- 
straining power (all of it in the case of z ~ 1) 
originates from k < 0.2 h Mpc -1 . At higher red- 
shifts the transition in these halo clustering mod- 
els scales to higher k as described in Smith et al. 
(2002, Figure 15). 

2.2. Back-of-the-envelope calculation 

A simple calculation demonstrates that to de- 
lineate the oscillations in the galaxy power spec- 
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trum accurately, we need to survey a cosmological 
volume greater than that of the local Sloan Digi- 
tal Sky Survey (which we take as a uniform cone 
of 10 4 deg 2 to z = 0.2, i.e. V S ioan = 2 x 10 8 h~ 3 
Mpc 3 ). There arc two sources of statistical error in 
a power spectrum measurement: sample variance 
(the number of independent wavelengths 27r/fc of 
a given fluctuation that can fit into the survey vol- 
ume V) and shot noise (the imperfect sampling of 
fluctuations by the finite number of galaxies N). 

The error due to sample variance on a power 
spectrum measurement, averaged over a radial bin 
in fc-spacc of width Afc, is 

( a P y_ (2tt)3 1 

(e.g. Peacock & West 1992). This expression is 
derived from the density of states and the volume 
of fc-space enclosed by the radial bin. The initial 
factor of 2 is due to the fact that the density field is 
real rather than complex, thus only half the modes 
are independent. (For a realistic survey geometry 
equation 2 is only an approximation, as neighbour- 
ing fc- modes are correlated). Figure 1 indicates 
that for a statistically significant measurement of 
the acoustic peak at k w 0.135 h Mpc - wc require 
ap/P < 0.02 as the fractional height of this peak 
is about 4 per cent. Thus taking Afc = 0.015 h 
Mpc- 1 (= k A /4) wc find that V > 1.8Vsioan- 

The number of galaxies needed to render shot 
noise negligible in this case is N ^> 1 jP «lx 10 5 
(assuming 6=1 and using the input model power 
spectrum to read off P x V « 5 x 10 3 h" 3 Mpc 3 
at k w 0.135/1 Mpc" 1 ). 

2.3. Simulations 

We constructed detailed simulations to deter- 
mine the scale of survey required to measure the 
oscillation wavescale k A with a given accuracy. 
When measuring a power spectrum from a real- 
istic survey, the quantity we derive is actually the 
power spectrum convolved with the survey window 
function in fc-space, W(k), the Fourier transform 
of the window function in real space. To detect 
baryonic oscillations it is important that W(fc) is 
compact (non-zero only for fc <C fc^), otherwise 
the smoothing effect of the convolution seriously 
reduces the amplitude of the oscillations. This re- 
stricts us to relatively simple survey geometries. 



(Any remaining convolution has no effect on the 
observed oscillation wavescale k A ). 

We considered two different potential surveys: 
the first directed at redshift range 0.5 < z < 1.3 
and the second probing the range 2.5 < z < 3.5. 
These ranges are reasonable possibilities for fu- 
ture large-scale surveys. The lower-redshift survey 
could target either luminous star-forming galax- 
ies via 3727A [Off] emission, or luminous ellipti- 
cal galaxies using absorption features such as the 
4000A continuum break. In cither case, z = 1.3 
is a reasonable limit where these spectral features 
move out of optical wavebands. For the higher red- 
shift survey, the Lyman Break selection technique 
operates efficiently for 2.5 < z < 3.5 as the 912A 
break shifts through the optical /7-band. Both 
these redshift regimes have already been explored 
by small surveys of hundreds of galaxies (e.g. Le 
Fevre et al. 1995, Steidel et al. 1998) and it is 
now becoming feasible to extend these efforts to 
far larger surveys. 

For each redshift range 0.5 < z < 1.3 and 
2.5 < z < 3.5 we defined a survey volume V 
by an angular patch on the sky of size 9x9. 
We varied the volume V by changing 9 such that 
< V < 6Vsioan in each case (0 max = 24.9° 
and 17.3° respectively). For the z ~ 1 survey 
we restricted the power spectrum measurements 
to the Fourier range k < 0.2 h Mpc -1 , which en- 
compasses the first two detectable acoustic peaks. 
Non-linearities and redshift-space distortions are 
likely to become important for k > 0.2 h Mpc -1 , 
diluting the baryonic features. Our estimate of 
the linear/non- linear transition scale k n \(z = 1) w 
0.2 h Mpc -1 is based on using the model power 
spectrum to calculate the dimcnsionless variance 
of mass fluctuations, er 2 (i?), inside a sphere of co- 
moving radius R equivalent to a fluctuation of 
wavelength 2ir/k. We defined equivalent scales 
by R = ir /2k, matching a half- wavelength (a full 
wavecrest) to the diameter 2R of the spheres. A 
conservative estimate of k n i(z = 0) is 0.1 /i -1 Mpc 
(Meiksin, White & Peacock 1999) which corre- 
sponds to (T 2 (-7r/2fc n i) = 0.35. At z = 1, this 
same value of a 2 is produced if fc n i w 0.2 h Mpc -1 
(see Figure 1). For the z ~ 3 survey we extended 
the measured power spectrum range to fc < 0.3 h 
Mpc -1 , which contains effectively all the visible 
acoustic peaks. The linear /non-linear transition 
scale at z — 3 is fc n i w 0.5 h Mpc -1 . 
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We generated many different Gaussian realiza- 
tions of N galaxies from the model P(k) within the 
survey volumes V , using Fast Fourier Transform 
techniques. For simplicity we assumed a constant 
selection function across the survey volume (i.e. 
a constant galaxy number density). A more real- 
istic flux-limited survey will detect galaxies with 
a varying redshift distribution N(z), but optimal 
weighting techniques for deriving P(k) exist in this 
case (Feldman, Kaiser & Peacock 1994). We also 
neglected the evolution of P(k) with redshift over 
the survey depth, fixing the model power spectrum 
at the mean redshift of the survey (i.e. z e s = 0.9 
or z c g = 3). This approximation has a negligible 
effect on the results as we are only interested in 
the relative positions of the peaks and troughs in 
P(k), which are independent of redshift. 

For the z ~ 1 survey we assumed a linear bias 
factor 6=1 for the galaxies. The potential z ~ 
3 survey will most likely be targetted at Lyman 
Break galaxies, which are known to be strongly 
biased tracers of mass (Steidel et al. 1998); in this 
case we assumed a linear bias parameter 6 = 3. As 
noted above, the main effect of this assumption is 
to reduce by a factor 6 2 the number of galaxies N 
required to render shot noise negligible. 

We measured the (noisy) power spectrum of 
each Gaussian realization using the standard 
method (see e.g. Hoyle et al. 2002). The mea- 
surements for the set of realizations allow the 
statistical error in each spatial frequency bin to 
be derived, and the ensemble incorporates any 
correlations that exist between adjacent bins. For 
each (N, V) simulation we generated 400 realiza- 
tions. For each realization we divided the mea- 
sured power spectrum by the smooth reference 
spectrum and fitted a decaying sinusoidal func- 
tion to the result, deducing a best-fit oscillation 
wavescale kA in A-space. We permitted the over- 
all amplitude A of the fitted sinusoidal function 
to vary but fixed the decay rate at an empirically- 
estimated value, so that we were simply fitting the 
two-parameter function: 



P{k) 



\+A k exp 



0.1/iMpc" 



1.4 



2irk 
kA 
(3) 




The power of 1.4 in the decay term originates 
from the Silk damping fitting formula presented 
in Eisenstein & Hu (1998, equation 21). Vary- 



Fig. 2. — Power spectrum measurement for a sim- 
ulated survey of N = 2 x 10 6 galaxies over a vol- 
ume V = 6 Vsioan at redshift z ~ 1. The power 
spectrum is divided by the zero-baryon model in 
the upper panel and the smooth reference spec- 
trum in the lower panel. The solid line is the 
input (unconvolved) model power spectrum (i.e. 
Figure 1). The dashed line in the lower panel is 
the best fit to the data of a simple decaying si- 
nusoidal function (equation 3). Points are plotted 
at intervals of Ak = 0.01 h Mpc -1 , which are ap- 
proximately uncorrelated (see Meiksin, White & 
Peacock 1999). 



ing the decay length as well as the amplitude was 
found not to have a significant effect on the fitted 
values of fc^. For small k there is a phase shift in 
the sinusoidal term (Eisenstein & Hu 1998, equa- 
tion 22) but this only affects the acoustic peak at 
k w 0.025 h Mpc -1 (see Figure 1) which is not 
measurable due to sample variance. 

For illustration, Figure 2 displays the power 
spectrum measurement and the best fit of equa- 
tion 3 for the first realization of the case z *~ 1, 
V = 6y S ioa„, N = 2x 10 6 . 

The distribution of fitted oscillation scales over 
the different realizations allowed us to describe 
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the accuracy of the measurement by a quantity 
Afc^/fcA, where AkA is half the range enclosed by 
the 16 th and 84 th percentiles of the distribution 
(i.e. Ua ± AkA defines the 68 per cent confidence 
region). Figures 3 and 4 plot contours of the accu- 
racy AkA/kA in the parameter space of V and N 
for the two potential surveys at z ~ 1 and z ^ 3. 
The contours running parallel to the V~ axis at 
small N indicate dominance by shot noise; the 
contours running parallel to the N— axis at large 
N represent dominance by sample variance. The 
diagonal dashed lines shown in Figures 3 and 4 
provide a reasonably good representation of the 
most efficient observational strategies. 

In order to measure the oscillation wavescale 
Ica with a precision of 2 per cent with negligible 
shot noise, we require either ~ 10 6 galaxies cov- 
ering ~ 400 deg 2 at z ~ 1, or ~ 5 x 10 5 objects 
over ~ 150 deg 2 at z ~ 3. The smaller number 
of galaxies needed at z ~ 3 is due to the larger 
bias factor b = 3 which boosts the amplitude of 
the density fluctuations, more than compensating 
for the smaller growth factor D\{z). In addition, 
a given angular field covers more volume at z ~ 3 
than at z ~ 1. However, Lyman Break galaxies 
will require longer exposure times. The bright- 
est z ~ 3 objects of Steidel et al. (1998) typi- 
cally required ~ 2 hours of integration on a 10m 
telescope to acquire a secure redshift for apparent 
magnitude R ~ 24. An L* galaxy at z ~ 1 has 
brightness R ~ 23 and should only require ~ 30 
minutes exposure. 

3. Measuring dark energy with baryonic 
oscillations 

In this section we investigate how accurately a 
detailed measurement of the baryonic oscillations 
can constrain the dark energy parameter w, which 
we assume does not vary with redshift. In order 
to measure P(k) from a galaxy redshift survey we 
must convert redshifts to co-moving co-ordinates 
(in h' 1 Mpc), assuming values for Q m and w. The 
expected wiggle wavescale fc^ is determined by the 
sound horizon before recombination; it is a func- 
tion of fi m , Qb and h (see equation 1 or the fit- 
ting formula of Eisenstein & Hu 1998, equation 
26). Incorrect cosmological parameters will lead 
to a distortion of the measured power spectrum 
so that the deduced value of fc^ is not consistent 




N / 1 

Fig. 3. — Fractional accuracy AkA/kA with which 
the wavescale of the baryonic oscillations in k- 
space can be measured at redshift z ~ 1, as a 
function of the number of galaxies N (as a frac- 
tion of 10 6 ) and the survey volume V (as a fraction 
of the Sloan volume Vsioan = 2 x 10 s h,- 3 Mpc 3 ). 
Contours are shown corresponding to (beginning 
in the bottom left-hand corner) AkA/kA — 10%, 
5%, 3% and 2%. The positions of the 2dF and 
Sloan surveys are marked on the plot for com- 
parison. This does not accurately represent their 
precision in measuring fc^, because the linear 
regime at redshift zero extends only to k w 0.1 h 
Mpc -1 (whereas the simulations assume k n \ = 
0.2 h Mpc" 1 ). Also, the Sloan (LRG) sample will 
possess a higher linear bias factor b ss 1.6 (whereas 
the simulations assume 6=1), which enhances the 
"effective" value of N by a factor b 2 . Hence the 
Sloan LRG sample may measure the position of 
the acoustic peak at k s» 0.075 h Mpc -1 to an ac- 
curacy of ~ 5% (Eisenstein, Hu & Tegmark 1998). 
The diagonal dashed line indicates the most effi- 
cient observational strategies: fewer galaxies will 
result in shot noise domination, and more galax- 
ies will be "wasted" . The dashed line corresponds 
to a surface density ~ 2400 galaxies deg -2 in this 
case. 



with the theoretical expectation. 

For the purposes of this initial study we as- 
sumed that the current values of the matter den- 
sity and Hubble's constant are known precisely, 
Q m = 0.3 and h = 0.7. We considered a range 
of values for the baryon density To create 

a concrete example, we assumed a fiducial value 
w = — 1 for the dark energy and investigated how 
accurately we could recover this value from the 
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Fig. 4. — The same plot as Figure 3 for simulated 
surveys at redshift z ~ 3. Note the different scale 
on the N— axis. Contours are shown correspond- 
ing to (beginning in the bottom left-hand corner) 
Ak A /k A = 10%, 5%, 3%, 2% and 1.5%. The diag- 
onal dashed line corresponds to a surface density 
» 3400 galaxies deg~ 2 . 

simulated surveys. 

3.1. Back-of-the-envelope calculation 

For a flat geometry, the fundamental relation 
between co-moving distance x and redshift z can 
be written in the form 

dx c 1 

dz H nU 2 ^(i + z f + (fi-i _ i)(i + z )3(i+w) 

(4) 

For values of w near —1, the second term inside 
the is small for z > 1 and the zeroth-order 

dependence is x oc i? f2 m . This conveniently 
cancels the zeroth-order dependence of the sound 
horizon scale on f2 m and Hq in equation 1. Thus 
our cosmological test has reduced sensitivity to un- 
certainties in f2 m and Hq- We quantify this in 
Section 4 below. 

Suppose fl m — 0.3 and the true dark energy 
parameter is Wtme = — 1- Consider a measuring 
rod located at z = 1. If the assumed cosmology 
is Wass = —0.9 (but constant with redshift), then 
using equation 4 the length distortion of the rod is 
dx' /dx — 0.975 if oriented radially (i.e. its ends at 
fixed redshifts) and x'/x = 0.980 if oriented tan- 
gentially (i.e. its ends at fixed angular positions). 
Thus the zeroth order effect is a re-scaling in the 
length of the rod by ~ 2% and the first order effect 
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Fig. 5. — The length distortion of a rod as a func- 
tion of redshift, supposing the true cosmology is 
f^m = 0.3, Wtme = — 1 and the assumed cosmol- 
ogy is Vt' m = 0.3, Wass = —0.9. The dashed and 
solid lines illustrate respectively the distortion if 
the rod is oriented radially (i.e. dx 1 /dx) and tan- 
gentially (i.e. x' d9/xd9 = x'/x). Thus it can be 
seen that the primary effect of assuming an incor- 
rect value of w is a re-scaling of distances away 
from their true values. 

is a radial/transverse shear. The full redshift de- 
pendence of these quantities is illustrated by Fig- 
ure 5. 

This distortion of spatial scales carries over di- 
rectly into fc-space, and suggests that a measure- 
ment of ftA with 2 per cent precision at z ~ 1 
translates into a measurement of w with an ac- 
curacy Ait; w 0.1. The corresponding figures at 
z = 3 are dx' /dx = 0.991 and x'/x = 0.981; the 
smaller distortion in dx suggests that the method 
is slightly less powerful when applied at z ~ 3. 

We note that in detail, when combined with 
CMB observations, the assumption of constant 
fl rn in this comparison is not quite right. The 
CMB power spectrum accurately constrains the 
angular diameter distance to the surface of last 
scattering, which removes one degree of freedom 
in the parameters: we cannot vary f2 m and a con- 
stant w completely independently. In addition, 
observations by the Planck satellite will accurately 
determine the quantity TO /i 2 within a few years. 
However, in this initial treatment, we consider our 
experiment as independent and suppose that f2 m 
is specified by external datasets. 

It is interesting to contrast the baryonic oscil- 
lations method presented here with the Alcock- 
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Paczynski test. The latter does not employ a stan- 
dard length-scale known in advance, but rather 
compares quantities parallel and perpendicular to 
the line-of-sight. Physically, it is sensitive to dis- 
tortions in dx/x. Numerically, this amounts to 
the difference between the solid and dashed lines 
in Figure 5 ( %M - 1 ~ % - ^ for small distor- 
tions). The baryon wiggles method uses a known 
ruler and is therefore independently sensitive to 
distortions in x and dx, through the parallel and 
perpendicular components of the power spectrum 
(although in this initial treatment we average the 
power spectrum over angles). Numerically, this 
is equivalent to the absolute deviation from unity 
of the solid and dashed lines in Figure 5. Thus 
the "lever arm" for detecting any departures from 
w = — 1 is a factor ~ 4 higher at z = 1 for our 
method and, in contrast to the Alcock-Paczynski 
test, it does not depend on small-scale non-linear 
clustering details. 

3.2. Simulations 

We simulated surveys in the redshift ranges 
0.5 < z < 1.3 and 2.5 < z < 3.5, as before, 
using the same measured ranges of k and linear 
bias factors. In each case we defined the survey 
volume by a 20° x 20° angular patch. This gen- 
erated a volume V = 3.9 Vsioan for the z ~ 1 sur- 
vey and V = 8.0 Vsioan at z ~ 3. We assumed 
redshift catalogues of N — 10 6 galaxies at z ~ 1 
and N = 5 x 10 5 galaxies at z ~ 3, such that 
the measurements were not limited by shot noise. 
These examples are somewhat arbitrarily chosen, 
but are indicative of future large-scale redshift sur- 
veys. We selected the galaxies from a redshift dis- 
tribution that populated the survey volume ap- 
proximately uniformly. 

We allowed the baryon density to vary over the 
range 0.1 < fib/fi m < 0.2 (centred on the re- 
sult flb/fl m = 0.15 assumed in Section 2). For 
each corresponding power spectrum we generated 
a set of galaxy realizations and converted these 
to redshift catalogues supposing a "true" cosmol- 
ogy w = — 1. We then re- measured P(k) for each 
realization for a range of "assumed" cosmologies, 
— 1.3 < w < —0.7. There are some theoretical 
reasons to suppose w > — 1, but we do not impose 
this restriction on our analysis (Caldwell 1999). 

For each assumed value of w we used 400 differ- 
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Fig. 6. — Power spectrum measurement for a sim- 
ulated survey with the same parameters as Fig- 
ure 2, except that the value w = —0.8 has been 
incorrectly assumed. The wavescale of the fitted 
function (the dashed line) is spuriously distorted 
compared to the theoretical expectation (the solid 
line). 

ent galaxy realizations to derive a distribution of 
fitted wavescales in fc-space. If the assumed value 
of w differs from the "true" cosmology w = — 1, 
then the fitted wavescales will be distorted (see 
Figure 6) and their distribution will not be centred 
on the expected wavescale (computed from Eisen- 
stein & Hu 1998, equation 26). This allowed us 
to assign a likelihood to each value of w, based on 
the position of the expected wavescale in the mea- 
sured distribution. For example, 68 per cent of 
the fitted wavescales lie between the 16 th and 84 th 
percentiles of the measured distribution. Thus if 
the expected wavescale lies at the 16 th percentile of 
the distribution, the value of w is rejected with "ler 
significance" . Likewise, if the expected wavescale 
lies at the 2.5 th percentile of the distribution, the 
value of w is rejected with "2er significance". 

Figures 7 and 8 display the final results for the 
z ~ 1 and z ~ 3 surveys: la and 2a confidence 
levels for the value of w as a function of fib- For 
the fiducial value of flb/Q m = 0.15 we can in- 
deed measure w to a la accuracy of ±0.1, bearing 
out our back-of-the-envelope calculation in Section 
3.1. This of course assumes the fiducial model 
w = —1; the error bar will have some dependence 
on the "true" value of w, which we do not investi- 
gate here. 

For comparison, we extended the angular area 
of the simulated z ~ 1 survey by a factor of five, 
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Fig. 7. — Constraints on the dark energy param- 
eter w produced by measuring the baryonic oscil- 
lations in the galaxy power spectrum, for a poten- 
tial large-scale redshift survey of 400 deg 2 at z ~ 1 . 
We assumed that the linear power spectrum can be 
successfully measured for k < 0.2 h Mpc" 1 . The 
survey parameters are V — 3.9Vsi oan , N = 10 6 , 
6 = 1. We assumed a "true" cosmology w = — 1 
and plot Oct, la and 2ct confidence limits on the 
measured value of w for varying SI;,. These con- 
tours correspond to the theoretical value of Ua 
lying at the 2.5 th , 16 th , 50 th , 84 th and 97.5 th 
percentiles of the observed wavescale distribution. 
The 50 th percentile is slightly displaced from w = 
— 1 because the full Eisenstein & Hu (1998) trans- 
fer function is used in the simulations and the 
approximate analytic formula (Eisenstein & Hu 
1998, equation 26) is used for the theoretical pre- 
diction. 

generating N = 5 x 10 6 galaxies over a volume 
V = 30Vsioan- In this case, the wiggle wavescale 
and dark energy parameter can be measured with 
respective accuracies Afc^/fc^ « 0.8 per cent and 
Aw w 0.04 (assuming f2j,/Q m = 0.15). 

4. Discussion 

Baryonic oscillations offer a measurement of 
the dark energy equation of state that comple- 
ments other methods. The precision obtainable 
from the surveys simulated above, Aw ~ 0.1 
(68% confidence), is comparable to that achieved 
by combining together all existing relevant ob- 
servations. The current high-redshift supernovae 
data requires w < —0.55 at 95% confidence for 
a flat geometry (Garnavich et al. 1998), and 
the latest (WMAP) CMB measurements, taken 
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w 

Fig. 8. — The same plot as Figure 7 for a simulated 
survey of 400 deg 2 at z ~ 3, supposing that the lin- 
ear power spectrum can be successfully measured 
for k < 0.3 h Mpc -1 . The survey parameters are 
^ = 8.0y sloan , iV = 5x 10 5 , 6 = 3. 

together with other datasets, imply w < —0.78 
at 95% confidence (Spergel et al. 2003). More- 
over, if deviations are to be detected from the 
w = — 1 paradigm, observations are required at 
z ~ 1 — * 3 where we expect the corresponding 
cosmological effects to be strongest: Universal ex- 
pansion switches from deceleration to acceleration 
at this epoch. 

A high-quality sample of 2000 high-redshift su- 
pernovae distributed from z = 0.2 to z = 1.7, 
provided by the SNAP satellite (Aldering et al. 
2002; http://snap.lbl.gov), promises to mea- 
sure a constant w to a precision Aw = 0.05 (68% 
confidence) and constrain its time variation. The 
latter is further improved when combined with 
the accurate measure of fl m h 2 provided by Planck 
satellite observations of the CMB (Frieman et al. 
2003). The baryon oscillations and supernovae 
methods of probing the dark energy should power- 
fully complement each other. Supernovae provide 
a "standard candle" approach (i.e. a measurement 
of the luminosity distance) , whereas baryon oscil- 
lations provide a "standard rod" approach, mea- 
suring this rod in both the transverse and radial 
directions. The transverse rod measures the angu- 
lar diameter distance. In conventional cosmology 
and physics, luminosity distance and angular di- 
ameter distance are both related to x(z) in equa- 
tion 4 by (1 + z) factors. However, this is not 
true of all scenarios; for an exotic alternative in- 
volving photons oscillating into undetectable light 
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axions see Csaki et al. (2002). The baryonic os- 
cillations additionally measure the standard rod 
radially, yielding a measurement of dx/dz which 
is equivalent to the Hubble parameter H(z). To 
the best of our knowledge the baryonic oscillations 
provide the only direct route to H(z) at high red- 
shift. 

The efficacy of the supernova approach is de- 
termined by how accurately systematic errors can 
be controlled. These systematic astrophysical ef- 
fects include the possible presence of non-standard 
host galaxy extinction, the evolution of supernovae 
with redshift, the effect of gravitational lensing, 
and variable Milky Way dust extinction. It is im- 
portant to correct the supernova magnitudes by a 
light-curve factor (Perlmutter et al. 1997) and/or 
colour/spectrum factor (Riess et al. 1998). In 
contrast, the acoustic oscillations method is not 
limited by systematic errors, but by statistical un- 
certainties due to sample size. 

The two approaches also differ in their re- 
quirement of a local calibration. In the super- 
nova method, cosmological parameters follow from 
the relative brightnesses of supernovae at differ- 
ent rcdshifts: hence the Nearby Supernova Fac- 
tory (http://snfactory.lbl.gov) plans a de- 
tailed study of several hundred low-redshift super- 
novae. In contrast, the measuring-rod length for 
the acoustic peaks is given directly by fundamen- 
tal physics (fl m h 2 and flbh 2 ) and does not require 
additional observations at z = 0. 

The use of counts of galaxies or clusters as 
probes of dark energy also depends on understand- 
ing critical systematics such as the variation with 
redshift of the intrinsic number density of the ob- 
jects, and the ability to identify systems of com- 
parable mass at different redshifts. 

The approach presented here is essentially a ge- 
ometrical method, similar in conception to the 
Alcock-Paczynski test. A classic application of 
the Alcock-Paczynski test is to small-scale galaxy 
clustering parallel and perpendicular to the line- 
of-sight (Ballinger et al. 1996). However, in this 
case it is very difficult to distinguish the geomet- 
rical distortions imprinted by an incorrect cosmol- 
ogy from redshift-space distortions. In contrast, 
baryonic oscillations probe clustering on much 
larger scales, where redshift-space distortions are 
insignificant. 



In our very empirical approach, the measure- 
ment of w depends only on knowledge of the po- 
sitions of the acoustic peaks and troughs in the 
power spectrum (and not on the detailed shape 
and amplitude of the power spectrum). It is 
the positions of these peaks, our "standard rod", 
which are the least model-dependent aspects of the 
power spectrum. Detailed models exist for physics 
before recombination, tested by increasingly ac- 
curate observations of the CMB power spectrum. 
The subsequent physics governing the growth of 
structures lies firmly in the well-understood linear 
regime for the scales of interest. 

However, three important issues have been ig- 
nored in this initial investigation and will be con- 
sidered in a future study (see also Seo & Eisen- 
stein, in prep.). The most important unexplored 
effect in the use of baryonic oscillations is the 
possible existence of more complicated galaxy bi- 
assing schemes (a non-linear bias that may evolve 
with redshift). This said, linear scale-independent 
bias should not be a bad approximation on large 
scales (e.g. Peacock & Dodds 1994). Secondly, 
we have restricted our attention to models with a 
constant, redshift-independent value of w. Our 
approach can be simply extended to constrain 
models of varying w, such as the common pa- 
rameterizations w(z) = wq + u>i z or w(z) = 
wo + w a z / (I + z) which will be the subject of a 
future study. Thirdly, we have neglected the un- 
certainty in the value of fi m , which clouds precise 
knowledge of the expected positions of the acous- 
tic peaks. 

It is easy to estimate the implications of uncer- 
tainty in Q m for our method, assuming a "true" 
cosmology fl m — 0.3, w = — 1. We saw in Sec- 
tion 3.1 that a shift Aw = 0.1 produces a radial 
distortion (in dx/dz) of w 2 per cent at z = 1. 
Using equation 4, a similar distortion is produced 
by a shift Af2 TO w 0.02. However, the sound hori- 
zon is also changed by w 1.5 per cent in the same 
sense (using the fitting formula of Eiscnstein & Hu 
1998, equation 26), thus the effects largely can- 
cel out. The "overall" radial distortion (relative 
to the sound horizon) becomes 2 per cent when 
Af2 m 0.07, which is reassuring: combining the 
WMAP CMB results with other astronomical data 
sets produces a current measurement of fl m with 
accuracy Af2 m w 0.04 (Spergel et al. 2003), and 
this situation should improve in the light of new 
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data. However, the z ~ 3 survey demands much 
tighter knowledge of Ct m than the z ~ 1 survey, 
because the (1 + z) 3 term inside the in equa- 
tion 4 starts heavily winning and the distortion 
due to a shift in w decreases. The result is that at 
z = 3, the shift Af2 m w 0.01 produces the same 
overall radial distortion as Aw = 0.1. Thus the 
application of this method at z = 3 requires more 
accurate knowledge of fi m than is achievable from 
current data. 

It is worth considering whether an imaging sur- 
vey using photometric redshifts could accurately 
measure the acoustic peaks (see also Eisenstein 
2002). Let us introduce a Gaussian scatter a z = 
(l + z)(To in each galaxy redshift. This corresponds 
to a smearing in co-moving co-ordinates of mag- 
nitude oy = a z (dx/dz). Taking a concrete exam- 
ple, if (T = 0.03 (a realistic estimate of the best 
possible current accuracy of photometric redshifts, 
e.g. Wolf et al. 2003) then oy w 100 hr 1 Mpc at 
z — 1. The effect on a spherically symmetric P(k) 
in (k x ,k y ,k z )-space is a damping of all modes by 
a factor exp(— k z a^) (in the flat-sky approxima- 
tion), where the z-axis is the line-of-sight. Con- 
sidering the modes in a spherical shell of radius 
k = 0.2 h Mpc -1 , the majority are too heavily 
damped to provide signal. The surviving useful 
modes (with \k z \ < 1/oy w 0.01 h Mpc -1 ) are lo- 
cated within a "ring" in /c-space, corresponding to 
a reduction in the number of useful modes by a 
factor ka r « 20. To achieve the same power spec- 
trum precision on this scale k, the sky area of the 
imaging survey (i.e. density of states in fc-space) 
must be increased by this same factor w 20 such 
that the ring contains the same number of modes 
as the original spherical shell. 

However, the most significant disadvantage of 
the photometric redshifts approach is one loses 
the capacity to decompose the power spectrum 
into its radial and tangential components, which 
separately constrain the Hubble constant and an- 
gular diameter distance through distortions in dx 
and x (Eisenstein 2002). In our initial "proof-of- 
concept" treatment, we have averaged the power 
spectrum over angles and not pursued these sep- 
arate constraints, which certainly merit more de- 
tailed calculation (see Seo & Eisenstein, in prep.). 
As noted above, the radial component of P(k) is 
damped by a factor exp (— fc 2 cr 2 ), which becomes 
unacceptably small unless ay < 1/fc or (taking 



k = 0.2/i Mpc -1 ) cr < 0.001 at z = 1. Spec- 
troscopic resolution R <~ 1000 is thus required for 
sufficiently accurate redshifts. 

Finally, how realistic is it to carry out such 
redshift surveys in practice? Our simulated sur- 
veys involve ~ 10 6 galaxies over ~ 400 deg 2 of 
sky. The exposure times for spectra of these ob- 
jects would be of order 1-2 hours on 8-metre class 
telescopes, based on existing data (e.g. Steidel et 
al. 1998). An 8-metre telescope with a 1° diam- 
eter field of view, capable of taking spectra of up 
to 3000 galaxies simultaneously, could accomplish 
such a survey in 60-120 nights. Such an instru- 
ment is eminently feasible and has been proposed 
several times (Glazebrook 2002). In particular, a 
detailed concept design has been put together for 
the Gemini telescopes (Barden 2002). The observ- 
ing time and data volume of such a survey is very 
similar to that of the 2dF Galaxy Redshift Survey 
(Colless et al. 2001). 

5. Conclusions 

This initial study has demonstrated that, under 
a simple set of assumptions, baryonic oscillations 
in the galaxy power spectrum may be used to mea- 
sure accurately the equation of state of the dark 
energy. For example, a survey of three times the 
Sloan volume at mean redshift z ~ 1 can measure 
the sound horizon (i.e. the wiggle scale) with an 
accuracy of 2 per cent, and the parameter w (as- 
sumed constant) with precision Aw = 0.1, assum- 
ing f2b/f2 m = 0.15. This probe of the dark energy 
is complementary to the supernova method, with 
an entirely different set of uncertainties dominated 
by sample size rather than systematic effects. It 
would provide the second independent pillar of ev- 
idence for the current epoch acceleration of the 
Universal expansion. The constraints on w be- 
come tighter as the survey volume is enlarged. We 
conclude that delineation of the baryonic peaks, 
and their use as a standard cosmological ruler to 
constrain cosmological parameters, should form an 
important scientific motivation for the next gener- 
ation of galaxy redshift surveys. 
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